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We analyze the effects of Zeeman degeneracies on the long-range interactions between like Rydberg 
atoms, with particular emphasis on applications to quantum information processing using van der 
Waals blockade. We present a general analysis of how degeneracies affect the primary error sources 
in blockade experiments, emphasizing that blockade errors are sensitive primarily to the weakest 
possible atom-atom interactions between the degenerate states, not the mean interaction strength. 
We present explicit calculations of the van der Waals potentials in the limit where the fine-structure 
interaction is large compared to the atom-atom interactions. The results are presented for all 
potential angular momentum channels invoving s, p, and d states. For most channels there are 
one or more combinations of Zeeman levels that have extremely small dipole-dipole interactions 
and are therefore poor candidates for effective blockade experiments. Channels with promising 
properties are identified and discussed. We also present numerical calculations of Rb and Cs dipole 
matrix elements and relevant energy levels using quantum defect theory, allowing for convenient 
quantitative estimates of the van der Waals interactions to be made for principal quantum numbers 
up to 100. Finally, we combine the blockade and van der Waals results to quantitatively analyze 
the angular distribution of the blockade shift and its consequence for angular momentum channels 
and geometries of particular interest for blockade experiments with Rb. 
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I. INTRODUCTION 

The interactions between ultracold Rydberg atoms 
have been of interest for some time now, beginning exper- 
imentally with studies of resonant energy transfer [If, Q ■ 
The strong, long-range interactions between Rydberg 
atoms are recognized as being extremely important for 
understanding these phenomena, such as the evolution of 
clouds of ultracold Rydberg atoms into ultracold plasmas 
and vice versa [3, 4J. Compelling theoretical concepts for 
exploiting Rydberg-Rydberg interactions for quantum 
information processing using single-atom qubitsj^] and 
atomic ensembles [6fl have stimulated further experimen- 

work. 
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One of the key ideas behind the potential for co- 
herent quantum information applications using Rydberg 
atoms is the concept of dipole blockade [6]. If an en- 
semble of N ground state atoms is resonantly driven by 
a laser tuned to a Rydberg state, the excitation pro- 
ceeds first through states with a single atom excited, 
then through states with two atoms excited, and so on. 
If the Rydberg-Rydberg interaction energies are suffi- 
ciently strong, excitation of the multiply-excited atom 
states will be greatly suppressed, and only a single atom 
will be excited at a time. Excitation of multiple atoms 
is blocked by the dipole-dipole interactions between the 
Rydberg atoms. In this blockade regime, an effective two- 
level system is realized between the ground state \ggg) 
(for N — 3) and the entangled symmetric excited state 

I s ) = (\99 e ) + \9 e 9) + l e <7S))/\/3- This two-level system 
has an effective light-atom coupling that is a factor of 
^/N larger than the light-single-atom coupling, and is 



promising for a wide variety of quantum manipulation 
applications. 

Some of the issues associated with realizing dipole 
blockade are illustrated in the energy level diagram for 
the case N = 3 in Figure [TJ The atoms are assumed to 
be placed symmetrically along a line so the dipole-dipole 
interactions between atom pairs are A12 = A23 > A13. 
Blockade is effective if and only if each of the dipole- 
dipole interaction strengths is much greater than the 
atom-light coupling. In that case, all the two- and three- 
atom excited states are out of resonance with the laser 
and an effective two-level system results. The influence of 
the two-atom states is limited to a small AC-Stark shift 
of the resonance frequency. However, if A13 is insuffi- 
ciently large then the doubly-excited state \ege) becomes 
appreciably populated. In addition, \ege) is also reso- 
nantly coupled to other singly-excited states that do not 
couple optically to the ground state. Thus the primary 
errors that enter the blockade process are the produc- 
tion of doubly-excited states [1] and singly-excited states 
outside the desired two-level system. 

For the case of non-degenerate Rydberg states, these 
errors have been discussed before [y, |24|. In Section HT1 
of this paper, we extend the analysis to the more real- 
istic case of degenerate alkali Rydberg states interacting 
via van der Waals forces. The types of errors are qualita- 
tively the same, but we show explicitly how to account for 
degeneracies in the Rydberg energy levels. We show that 
the figure of merit for blockade is primarily determined 
by the weakest of the many potential curves that emanate 
from a given degenerate Rydberg state. In many cases 
there are one or more of these degenerate states with 
nearly zero van der Waals interaction and hence weak 
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FIG. 1: (color online) Energy levels of 3 atoms arranged in a 
line, with dipole-dipole interactions A12 = A23 > A13. The 
dashed lines are the energy level positions if the dipole-dipole 
interactions were zero, and the arrows denote the successive 
excitation of the multiple atom states by the resonantly-tuned 
laser. Solid lines indicate pairs of states that are coupled 
by allowed electric dipole transitions. All these transitions 
are off-resonant save the \ggg) — \s) pair that consititute an 
effective two-level system in the presence of blockade. 



blockade, despite the state-averaged van der Waals inter- 
action being substantial. These "Forster-zero" states [13] 
can be thought of as resulting classically from precession 
of the atomic angular momenta due to the dipole-dipole 
interaction. These states did not appear in Ref. due 
to the use of non-degenerate perturbation theory, effec- 
tively ignoring the ability of the dipole-dipole interaction 
to change the orientation of the individual atomic angu- 
lar momenta. 

With an eye toward applying these ideas to Rydberg 
states in the 30 < n < 100 range of interest for quantum 
manipulation of atoms at optically resolvable distances, 
we are led in Section IIIII to consider van der Waals in- 
teractions at distances where they are weak compared to 
the fine-structure splitting, a limit rarely treated in the 
literature [13, [H, [30] - We calculate the van der Waals 
eigenvalues for 21 different angular momentum channels 
likely to be encountered for low-angular momentum Ry- 
dberg states, and present the numerical values of ma- 
trix elements and energy defects so that quantitative es- 
timates of blockade can be made by others. We follow 
this in Section IIVI with two concrete calculations of the 
blockade frequency shifts to be expected for the case of a 
linear geometry, and the angular distributions for inter- 
esting cases. 

We note that in this paper we are using the term 
"blockade" in a more restrictive sense than has been used 
in a number of experiments to date. It has become com- 
mon to refer to a suppression or reduction in the number 
of excited Rydberg atoms due to atom-atom interactions 



as blockade. However, such effects do not necessarily im- 
ply the type of blockade necessary to produce collective 
singly-excited quantum states that are of primary inter- 
est for this paper. 



II. ERRORS IN QUANTUM MANIPULATIONS 
USING BLOCKADE 



In our previous work [24f | on single atom and single 
photon sources, we assumed spatially uniform Rabi fre- 
quencies and a spherical distribution of atoms. We also, 
in order to keep things simple, assumed the effective Ry- 
dberg shift for each atom was the same, and assumed a 
non-degenerate Rydberg state. We wish to remove these 
assumptions in order to more accurately describe situa- 
tions that will be faced in Rydberg blockade experiments. 

We will assume that the blockade shifts are large 
enough to limit a cloud of N atoms to at most 2 excited 
atoms. Thus the atomic cloud can be in the possible 
states 



\g), I'yk), and \ipkl) 



(1) 



representing respectively all the atoms in the ground 
state, the fcth atom in the singly-excited Rydberg state 7, 
and the kth and Ith atoms in the doubly-excited Rydberg 
state ip. 

We assume that the coupling of the excitation light to 
the atoms can be represented by an effective excitation 
operator 



(2) 



that takes atoms at various positions k from the ground 
state to a Zeeman sublevel \-fk) of the Rydberg state, 
and it also takes states with one Rydberg atom excited at 
position A: to a doubly-excited state. The Rabi coupling 
to the state \ jk) is 



(3) 



Here and in the remainder of the paper we use units with 

n = i. 



Generalizing Saffman and Walker [2 4j. the wavefunc- 
tion for the TV-atom ensemble is 

= Cg\g) + ^ C ; fe | 7 /c) + c*ki\<pkl) (4) 

rk tp,k<l 

where the doubly-excited states \ipkl) are eigenstates of 
the effective Rydberg-Rydberg Hamiltonian: 



H e s\ipkl) = A vk i\tpkl) 



(5) 



The index ip represents the various possible Rydberg- 
Rydberg shifts A^j-i for atoms k and I in the sample. 
For the specific situations treated in this paper, H e ff rep- 
resents van der Waals interactions and so rotational in- 
variance guarantees that A v ki depends only on the index 
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if and the distance R between atoms k and /. Likewise, 
the states \ipkl) for different atom pairs are independent 
of k and I when expressed in a coordinate system oriented 
along the interatomic axis; in the laboratory frame they 
are related to each other by rotations. We will not use 
these properties in the remainder of this section, so the 
results also apply to cases where the effective Rydberg- 
Rydberg interactions are modified by external laboratory 
fields that do not preserve their rotational invariance. 

It is convenient to introduce the wavefunction for a 
symmetric singly-excited state as follows: 
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and orthogonalize the remaining singly-excited states as 



|(7fc)x>=a 7fc (| 7 fe)-<s| 7 fc)|s)) 



(7) 



where a 7 fe is a normalization factor. These orthogonal- 
ized states do not couple directly to \g) via the light, 
but are populated only by coupling through the doubly- 
excited states. The collective Rabi frequency is defined 
to be 



ond excited states are 
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(11) 



(12) 



(13) 



where we have omitted the possible excitation of the 
doubly-excited states from the orthogonalized singly- 
excited states. 

We solve these equations as a successive approximation 
in the ratios fi/A. The doubly-excited amplitudes are of 
order £1/ A so to zeroth approximation we take c v ki — 
(perfect blockade approximation) and get 
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(14) 



which are the equations for standard Rabi flopping at col- 
lective Rabi frequency Qn- Explicitly, a resonant Rabi 
pulse of duration T beginning with c g = 1, has the solu- 
tion 
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where Slo is the rms single-atom Rabi frequency averaged 
over the Rabi frequencies of the individual atoms. 

The symmetric singly-excited state is coupled by the 
light to the ground state and the doubly-excited states. 
The matrix elements are 



(s\A^\g) 
(<pkl\A^\s) 
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where the factor of 2 disappeared because both k' = 
k,l' = I and k' = 1,1' = k terms contribute to the sum. 
The dimensionless overlap factor k^m gives the relative 
amplitude for exciting a particular pair kl of atoms to 
the doubly-excited Rydberg state (p. It depends on the 
experimental geometry (laser polarization, spatial varia- 
tion of intensity, etc.) and in particular depends on the 
relative orientation of the atom pair kl and the light po- 
larization. 

Using the above definitions, the Schrodinger equations 
for the ground, symmetric first-excited state, and the sec- 



cos(J1jvT/2) 
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In the approximation of perfect blockade, the system un- 
dergoes collective Rabi flopping without dephasing or loss 
of population to other states. 

Let us now consider the effectiveness of Rydberg block- 
ade. We start by calculating the probability of excitation 
of more than one atom. In the limit of large but finite 
Rydberg- Rydberg shift, we can make an adiabatic ap- 
proximation to Eq. (fl"3"|) to get 
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The probability of double excitation is 
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(17) 



(18) 



It is critical to note that, given relatively even excita- 
tion of the two-atom Rydberg states, it is an average of 
1/A^ fc ; that determines the blockade effectiveness. This 
means that Rydberg-Rydberg states with small van der 
Waals shifts are much more strongly weighted than those 
with large energy shifts. Let us define a mean blockade 
shift B via 



1 

B 2 



2 ^ K v> kl 
N(N — 1) ^ A 2 ,,., 

v ' ip-k<l fkl 



(19) 
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Then the probability of double excitation becomes 



Pi 



2MB 2 



(20) 



This shows that, for fixed f2jv, the probability of dou- 
ble excitation is virtually independent of the number of 
atoms in the ensemble. This does not contradict Fig. 3 
of Ref. 24j , where the plot assumed a fixed value of the 
single-atom Rabi frequency, not f2jy. We will evaluate B 
for cases of experimental interest below. 

It is is important to keep in mind that the blockade 
shift B depends on the polarization of the excitation light 
as well as the Zeeman structure of the state \g) through 
the overlap factor n v ki- We do not explicitly indicate 
these dependences in order to avoid a proliferation of 
subscripts. Explicit examples will be given in Section HVl 

In addition to the production of population of doubly- 
excited states, finite blockade also causes a frequency 
shift of the effective two-level system through virtual ex- 
citation of the doubly-excited states. Using the adiabatic 
approximation results for c v ki modifies Eq. (|12p to 
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The second term represents a shift in the resonance fre- 
quency of the effective two-level system. Defining a fre- 
quency shift factor 
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the resonance frequency of the two-level system is shifted 
by 
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At this level of approximation, the final blockade er- 
ror that can result is to transfer population out of the 
"computational basis" of \g) and |s) and into the other 
singly-excited states |(7fc)_i-)- The amplitude for these 
states obeys 
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in the adiabatic approximation (| 1 T[) . This does not sim- 
plify as nicely as Eq. (fl"9]l . but a good estimate of the 
probability of excitation of the singly-excited states out- 
side the computational basis is 



P[~P 2 



(26) 



Thus in estimating the blockade errors one should rough- 
ly double the estimate (|2"0"|) obtained from the blockade 
shift and the collectively enhanced Rabi frequency. 



III. VAN DER WAALS INTERACTIONS 
BETWEEN DEGENERATE RYDBERG ATOMS 

A. General Discussion 

The electrostatic dipole-dipole interaction between two 
non-overlapping atoms A and B that lie a distance R 
apart along the z-axis is 
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where a is the position of the electron on atom A and b is 
the position of the electron on atom B, and p = —p. The 
odd parity of the dipole operators a and b results in an 
initial 2-atom state \n aI aJ Am AnB^BjBms) being mixed 
with states of I a ± 1 and Jb ± 1. The individual total 
angular momenta ja and js may also change, consistent 
with dipole selection rules, as do the principal quantum 
numbers ua and ub- The total projection of the angular 
momentum along the z-axis M = rriA +mb is conserved, 
but the individual quantum numbers change by ±1 or 0. 

For the simplest version of blockade physics, we are 
interested in the case where the two atoms are being ex- 
cited to the same energy level, so that ua = ns = n, 
I A — Ib = I, 3 A = Jb =3- Then the dipole-dipole inter- 
action causes the reaction 



nlj + nlj — > n s l s j s + n t l f jt 



(28) 



with an energy difference between the final (s, t) and ini- 
tial two-atom states 



S = E(nJ s j s ) + E(n t l t jt) - 2E(nlj) 



(29) 



that we will call a Forster defect. At the largest dis- 
tances, a non-zero Forster defect causes the dominant 
long-range interaction between the atoms to be of the 
R~ 6 van der Waals type that arises from Vdd in second 
order. (The R~° quadrupole-quadrupole interaction is 
normally much smaller, as estimated in Appendix [A"l ) 
For atoms that are closer, the van der Waals interac- 
tion becomes large enough to mix the fine-structure lev- 
els together, particularly for the d levels. This occurs 
(for 30-80d levels) in the 0.8-8 /xm range of distances of 
interest for interactions between optically resolvable Ry- 
dberg atoms. For atoms at somewhat smaller distances, 
typically 0.5-5 /xm, the dipole-dipole interaction is com- 
parable to the energy differences between nearby states, 
so the interactions become resonant and vary as R~ 3 . 

There are some situations where the Forster defects 
are smaller than the fine-structure splitting, a notable 
example [TtJ being at 43d in Rb, where the reaction 
43d 5/2 + 43d 5/2 -> 45p 3/2 + 41/ has less than 10 MHz 
Forster defect while the 43d fine-structure splitting is 
150 MHz. In this case the transition to resonant dipole- 
dipole coupling occurs at longer range than fine-structure 
mixing. 
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As discussed above, the possibility of using Rydberg 
blockade for mesoscopic quantum manipulation depends 
critically on the most weakly interacting atoms, those 
that are furthest apart. It is therefore extremely impor- 
tant to understand the interactions of Rydberg atoms in 
the limit where the van der Waals interaction has not 
mixed the fine-structure. For most of the rest of this 
paper we will restrict our discussions to this limit. 

For a non-zero Forster defect, the energy shifts of the 
initial states are determined at long range by the effec- 
tive second-order perturbation operator or van der Waals 
interaction 

„ v- Vdd\st) (st\V dd 

-HvdW = > J (oUJ 

— Ost 

st " 

This operator is understood to act within the degenerate 
set of Zeeman sub-levels of the two-atom initial state. In 
general, the sum is over the various intermediate two- 
atom energy levels \st) that obey the selection rules of 
Vdd discussed above, and S st is the Forster defect for each 
channel with respect to the initial state. The indices s 
and t denote the full set of quantum numbers that spec- 
ify the intermediate states. In practice, this sum can be 
greatly simplified by noting that in most cases only states 
close in energy to the initial state have significant r ma- 
trix elements and therefore dominate the matrix elements 
of Vdd; and that typically one or two of these states have 
the smallest Forster defects. Thus there will usually be 
only a couple of intermediate states that give by far the 
biggest contributions to the van der Waals interactions. 

It is particularly important, as we will see, to properly 
account for the Zeeman degeneracy of the initial and in- 
termediate states. This degeneracy has the consequence 
that one cannot calculate the energy shifts by simply 
taking the expectation values of the van der Waals oper- 
ator of Eq. ([30]) . The dipole-dipole interaction couples an 
initial state with magnetic quantum numbers ttia, tub to 
intermediate states tua +p, tub —p which can then couple 
to a different Zeeman combination mA+p + q,mB — p — q 
where p and q range from —1 to 1. Thus the van der 
Waals interaction, being second order in Vdd, changes 
the magnetic quantum numbers of the individual atoms 
by up to ±2 units. It is therefore necessary to use degen- 
erate second-order perturbation theory [31| (Section 31) 
to understand the van der Waals interactions of Rydberg 
atoms. The diagonalization of the effective Hamiltonian 
of Eq. ([30]) leads to a range of eigenvalues, meaning that 
the strength of the van der Waals interactions depend 
strongly on the Zeeman sublevels. 

As we have shown previously [27j . many states which 
are nearly Forster resonant have linear combinations of 
Zeeman sublevels with zero dipole-dipole coupling. The 
existence of these states can be understood from angular 
momentum arguments. If the initial state atoms have 
individual angular momenta j, these can be coupled to- 
gether to make a total of 2j + 1 possible states of an- 
gular momentum J with M = 0. (Most of the Forster- 
zero states have M = so we restrict our argument to 



that case.) Similarly, if the dipole-dipole coupled states 
are constructed from angular momenta j s and jt, with 
js < jt, there are 2j s + 1 M = intermediate states. 
The condition for a "Forster zero" state ipp is 

2j 

(st|Vdd|V>F> = ( st \Vdd\ J) (J\*Pf) = (31) 

which must hold for all of the 2j s + 1 coupled states (st\. 
The Forster-zero condition is therefore a set of 2j s + 1 
equations in 2j + 1 unknowns. A Forster-zero solution 
exists for j s < j, and does not exist for j s > j. For j s = j 
there are not always precise zeros but in the cases studied 
here there is always a state with very small dipole-dipole 
coupling. 

The above arguments strictly hold for cases where a 
single channel dominates the dipole-dipole interactions. 
Additional channels can give some dipole-dipole cou- 
pling to the Forster zero states, but the states are still 
very weakly interacting, leading to potential problems for 
blockade applications. 



B. van der Waals Interactions with Fine Structure 

We now proceed to calculation of the dipole-dipole 
interaction at such long interatomic distances that the 
atomic fine-structure is not affected by the dipole-dipole 
interaction. We assume that all the intermediate states 
that are coupled to the initial states have the same angu- 
lar momentum structure, i.e. that a single virtual process 

nlj + nlj -> nj s j s + n t l t jt (32) 

occurs where various values of n s and n t may con- 
tribute but only a single value of l s ,js,h>jt- We de- 
note the initial states of the two atoms \nljmAnljrriB) 
in i in :> and the dipole-dipole coupled intermediate 
states \n s l s j s m s n t l t j t m t } = \m s m t ). 

The dipole matrix elements in Vdd can be written in 
terms of radial matrix elements and angular momen- 
tum factors using the Wigner-Eckart theorem [32| (Eq. 
13.1.5(40)): 

(n'l'j'm'\r p \nljm) = c££ 

x{j/ \ j,} {n>l>\\r\\nl) (33) 

where the reduced matrix element is 

(n'l'\\r\\nl) - v^TlC^o^'' = V^TTC^ 

x J rP n , v {r)P nl {r)dr (34) 

The radial wavefunctions P n i (r) can be calculated numer- 
ically using quantum defect theory or model potentials. 
It is convenient to define an operator M. that includes 
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all the angular momentum properties of Vdd for the in- 
termediate states of angular momenta j 8 , j t : 



(m s m t 
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(The matrix elements of M. are much simpler in a cou- 
pled basis (Appendix |5J rather than the product basis 
we are using here, but for connection to blockade physics 
the product basis is more convenient.) This allows the 
degenerate Hamiltonian due to the van der Waals inter- 
action to be written 

^6 . a>!\ \ i i . * C*6, 
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R 6 



V (36) 



where Cq depends only on the atomic energy level struc- 
ture and radial matrix elements: 
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(37) 



The operator T> — J\4^A4 contains all the angular mo- 
mentum properties of the states. Its (2j + l) 2 eigenvalues 
D v , when multiplied by Cq, give the long-range energies 
of the two-atom eigenstates: 



vdW 



\<P) 



R 6 



DM 



(38) 



The eigenvalues obey < D v < 1. The sign of the 
Forster defects determines the sign of Cg . In the case that 
channels of different angular momentum structure con- 
tribute significantly to the long-range interactions, the 
i/vdw matrices for each channel should be computed sep- 
arately, added together, and then diagonalized. 

We have calculated the eigenvalues (Table [J) and eigen- 
vectors [331 ] of T> for initial s, p, and d states with fine 
structure, corresponding to 23 different angular momen- 
tum channels. We will now discuss a few interesting cases 
before proceeding to using these results for blockade es- 
timates. 

We note that of the 23 different channels, 9 of them 
have at least 1 Forster zero state with zero dipole-dipole 
interaction. Another 7 have at least one state with 
an eigenvalue less than 0.05. All of these states are 
extremely weakly coupled by dipole-dipole interactions, 
and therefore are of limited use for blockade experiments. 

The remaining seven channels 
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(39) 



have minimum eigenvalues of 0.18 or greater. All of them 
have the property that the state being coupled to by the 
dipole-dipole interaction has larger angular momentum 
than the initial state. This is consistent with the argu- 
ment given above ( Section IlII Ap and in Ref. [27j that the 
preferred channels for blockade experiments have domi- 
nant intermediate channels where both atoms have total 
angular momentum j + 1. 

The three channels with no fine structure specified in 
the final state deserve special mention. Each of them 
have two allowed intermediate angular momentum chan- 
nels that contribute to the van der Waals interactions. 
When there is no near Forster resonance for one of those 
channels, to a good approximation the fine-structure in 
the final state can be neglected. An example of this will 
be given in Section [TV] 

The simplest case with zeros is P1/2+P1/2 ~^ s i/2+ s i/2- 
In this case the M — portion of the T> matrix is 



V 




(40) 



which has eigenvalues of and 16/81. The zero eigen- 
vector, — ("j"^)' nas zer0 contribution to the van- 
der- Waals interaction from the s + s states. 

For the case of p 3 / 2 +P3/2 ~> s i/2 + s i/2i the M = ±2 
and M = ±3 states have zero dipole-dipole coupling since 
the s + s states are limited to M = 0, ±1. Only 4 of the 
16 possible two-atom states have non-zero contributions 
to the van-der- Waals interaction from the s + s channel. 

An interesting case occurs in Rb due to the near res- 
onance of 43g?5/2 + 43g?5/2 — > 45^3/2 + 41/. The fine 
structure splitting of the f-states is small, so both the 
/ 5 / 2 and /V/2 states must be taken into account. In the 



approximation of zero / state fine structure splitting, we 
find that the eigenvalues of T> range from 0.836 down to 
0.0024, a factor of 350. The distribution of eigenvalues is 
shown in Figure [H 

Although the figure of merit for blockade is primarily 
determined by the energy shifts A v ki , the angular distri- 
bution and polarization dependence can play an impor- 
tant role for specialized geometries. These effects show 
up in the overlap factor 



(41) 



where \7k-fl} is the doubly-excited state that would be 
generated by the light in the absence of dipole-dipole in- 
teractions. The wavefunctions \ipkl) are simplest when 
expressed in a coordinate system aligned with the inter- 
atomic separation, while \jk"fl) is naturally represented 
in a fixed coordinate system. It is therefore convenient to 
rotate the van der Waals eigenstates to the fixed frame, 
denoted by primes, to get 



(ipkl\jk"fl) = 



(ipkUmkmAd 3 , d , 



X (™4 TO ; \-yk-yl) 



(42) 
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TABLE I: Relative interaction strengths for van der Waals interactions of Rydberg atoms, for various collision channels. The 
potential energy at distance R is the product CqD v /R 6 , which contains the effects of Zeeman degeneracy, with the overall Ce 
coefficient (Eq. I37[) for a particular channel that depends only on the energy level structure and radial matrix elements. Cases 
where the j quantum number is not included in the channel description are the sum over fine-structure components of the final 
state. 



Channel 




{^} 


Channel 


\M 


1 {D v } 


Sl/2 + 81/2 — » 


1 


{1.33} 


Sl/2 + S1/2 — > 


1 


{0.0988} 


P+P 





{1.33, 1.33} 


Pl/2 +P1/2 





{0.395, 0} 


Sl/2 + Sl/2 — » 


1 


{0.346} 


Sl/2 + Sl/2 — > 


1 


{0.543} 


Pl/2 +P3/2 





{0.444, 0.0494} 


P3/2 +P3/2 





{0.84, 0.444} 


Pl/2 +P1/2 — > 


1 


{0.0988} 


Pl/2 +P1/2 — » 


1 


{0.346} 


s l/2 + s l/2 





{0.395,0} 


s l/2 + G?3/2 





{0.444, 0.0494} 










3 


{0} 


Pl/2 +P1/2 — » 


1 


{0.543} 


P3/2 +P3/2 -» 


2 


{0,0} 


G?3/2 + ^3/2 





{0.84, 0.444} 


Sl/2 + Sl/2 


1 



{0.543, 0, 0} 
{0.84,0.444,0,0} 




3 


{0} 




3 


{0.267} 


P3/2 +P3/2 ~* 


2 


{0.08,0.00889} 


P3/2 +P3/2 -> 


2 


{0.48, 0.0533} 


S1/2 + d-j/2 


1 


{0.0622,0.0491,0.00322} 


Sl/2 + ^5/2 


1 


{0.64,0.0533,0.16} 







{0.0494,0.0178,0,0} 







{0.693,0.267,0,0} 




3 


{0.0128} 




3 


{0.0725} 


PS/2 + P3/2 -» 


2 


{0.00569, 0} 


P3/2 +P3/2 -» 


2 


{0.0672, 0.0587} 


G?3/2 + ^3/2 


1 


{0.00626, 0.00291, 0.00142} 


G?3/2 + ^5/2 


1 


{0.0654, 0.0608, 0.0189} 







{0.0178, 0.0149, 0.00217, 0.00182} 







{0.0675, 0.0657, 0.0328, 0.000416} 




3 


{0.269} 




3 


{0.354} 


P3/2 + P3/2 -> 


2 


{0.576,0.269} 


P3/2 +P3/2 -> 


2 


{0.635, 0.342} 


ds/2 + d 5 /2 


1 


{0.831,0.499,0.264} 


d + d 


1 


{0.857,0.561,0.33} 







{0.935,0.649,0.428,0.253} 







{0.948,0.699,0.501,0.321} 




3 


{0} 




3 


{0} 


G?3/2 + d 3 / 2 — » 


2 


{0,0} 


G?3/2 + ^3/2 — » 


2 


{0.08, 0.00889} 


Pl/2 +P1/2 


1 


{0.543, 0,0} 


Pl/2 +P3/2 


1 


{0.0622, 0.0491, 0.00322} 







{0.84,0.444,0,0} 







{0.0494, 0.0178,0,0} 




3 


{0.0128} 




3 


{0.267} 


^3/2 + G?3/ 2 — > 


2 


{0.00569, 0} 


G?3/2 + ^3/2 — > 


2 


{0.48, 0.0533} 


P3/2 +P3/2 


1 


{0.00626, 0.00291, 0.00142} 


Pl/2 + /5/2 


1 


{0.64,0.16,0.0533} 







{0.0178, 0.0149, 0.00217, 0.00182} 







{0.693,0.267,0,0} 




3 


{0.0725} 




3 


{0.269} 


^3/2 + G?3/2 — > 


2 


{0.0672, 0.0587} 


^3/2 + ^3/2 — > 


2 


{0.576,0.269} 


P3/2 + /5/2 


1 


{0.0654, 0.0608, 0.0189} 


/s/2 + /s/2 


1 


{0.831,0.499,0.264} 







{0.0675, 0.0657, 0.0328, 0.000416} 







{0.935, 0.649, 0.428, 0.253} 




5 


{0} 




5 


{0.343} 




4 


{0,0} 




4 


{0.503,0.137} 


C?5/2 + ds/2 — > 


3 


{0.269,0,0} 


C?5/2 + d 5 /2 — > 


3 


{0.642,0.263,0.0702} 


P3/2 +P3/2 


2 


{0.576,0.269,0,0} 


P3/2 + / 


2 


{0.747,0.373,0.145,0.0659} 




1 


{0.831,0.499,0.264,0,0} 




1 


{0.814, 0.442, 0.229, 0.0764, 0.0594} 







{0.935, 0.649, 0.428, 0.253, 0, 0} 







{0.836, 0.466, 0.243, 0.152, 0.00304, 0.00239} 


^5/2 + d 5 /2 — » 
/ + / 


5 
4 


{0.185} 2 {0.809,0.561,0.353,0.206} 
{0.397,0.193} 1 {0.938,0.703,0.505,0.332,0.208} 


3 


{0.619,0.375,0.201} {0.985,0.763,0.597,0.458,0.317,0.207} 
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FIG. 2: (color online) Range of van-der-Waals coefficients for 



the channel d 5 / 2 + ds/2 
lying nearly on the D v - 



-> P3/2 + j "5/2,7/2 • The M = point 
axis experiences little blockade. 



where the d's are Wigner rotation matrices evaluated at 
angle 9ki , the angle between the interatomic axis and the 
z-axis of the fixed coordinate system. 



C. Connection between van der Waals and Forster 
Regimes 

In the case that a single channel dominates, the eigen- 
states I ip) and eigenvalues D v of Ai^Ai can be used to 
analytically find the energies and eigenstates in the R~ 3 
Forster regime as well. The transition between van der 
Waals and Forster interactions occurs at a characteristic 
length scale of R c — (4C|/<5 2 ) 1 / 6 . A related analysis in a 
different context was given in Ref. [34| . 

The M. operator acting on a Forster eigenstate \<p) 
produces a unique vector \x v ) that is a superposition of 
the Zeeman sublevels of the coupled state: 

M\<p) = ^D~M (43) 

Operating on the left side with M. JA' we get 

MM*M\<p) = MD v \tp) (44) 
MM^\ Xv ) = DM (45) 

so \xip) is an eigenvector of M.M' with eigenvalue D v . It 



also follows that M^\xip) = 



,\ip), and (x v \x v ) 



1. 

Therefore the states \<p) and \x v ) form a closed two- level 
system under the influence of the dipole-dipole interac- 
tion. 

It is now straightforward to find the eigenstates and 
eigenvalues in the Forster regime as well. The Hamilto- 
nian matrix for the effective two-level system is 




(46) 



> 



4- 

2-. 


-2 H 

-4 

-6 




+43d c 



|M| = 
4 

3 
2 
1 





FIG. 3: (color online) The potential curves for 
43d 5/ / 2 + 43d 5/2 — > 45p 3/2 + 41/ in Rb, using S = -7.4 
MHz and C 3 = 1.98 GHz ^m 3 which give R c = 8.1 /Urn. At 
small R additional channels contribute so these curves are 
not accurate there. 



where 5 — E v ~ E x \s the Forster defect and Cq = Cf / S. 
The eigenvalues are 



V± (R) = -±^/52+4ClD v /RS 



and the eigenvectors are 



= cos9\(p) - sin^lx^) 
ip + = sm9\ip) +cos9\x v >) 



(47) 



(48) 



where tan 29 = -2C 3y ^/(5R 3 ). 

These considerations show that the long-range poten- 
tials for a + a — * b + b and the reversed b + b — ► a + a 
are anti-symmetric in energy about 6/2. If jb < j a , the 
channels with \M\ > 2jb have no dipole-dipole interac- 
tion. 

In the limit of strong dipole-dipole coupling, we have 



V± 



(49) 



The potential curves for 43d 5 / 2 + 43d 5 / 2 — > 45p 3 / 2 + 
41/, generated from Table fl] and the analytical formula 
above, are shown in Figure [3] 



D. Evaluation of van der Waals Interactions for Rb 
and Cs Rydberg States 

We now proceed to quantitatively evaluate the van der 
Waals interactions of Rydberg states that can be reached 
by one or two photon excitation from the ground state of 
neutral alkali atoms, with the restriction that we include 
only cases where both atoms are initially excited to the 
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same level. For each choice of excited state we give nu- 
merical values for the energy defects and the interaction 
strength for n = 70 as well as for values of n where reso- 
nances occur. Results are given for the two heaviest alkali 
atoms Rb and Cs. Due to the large hypcrfine splittings 
of the upperstates of the Dl and D2 lines in these atoms, 
they are the most promising candidates among the alkali 
metal atoms for quantum logic experiments which rely 
on a well resolved excited state hyperfine structure for 
qubit initialization and readout. 

Before discussing the cases individually we recall that 
the long range interaction strength of a particular chan- 
nel scales proportional to Cq = (i?^ s i?^*) /\S\ where 
we have introduced a shorthand notation 7 = {nlj} for 
the quantum numbers specifying the initial laser excited 
(7i) and Forster coupled (7 8) 7t) states. We have calcu- 
lated these matrix elements in several ways: using model 
potentials [35[, with quantum defect wavefunctions [36| . 
and using a semiclassical analytical formula[37]]. Results 
from the quantum defect wavefunctions and semiclassi- 
cal calculations typically agree to better than 1%, while 
the model wavefunction calculations differ from these by 
up to 10%. The numerical values in what follows were 
obtained from the quantum defect wavefunction calcula- 
tions. The radial integrals depend explicitly on n, I but 
also have an implicit dependence on j due to the de- 
pendence of the quantum defects on the fine structure 
level. This leads to a j dependence of as much as 10% 
in some cases. We report numerical values for the radial 
integrals corresponding to the particular fine structure 
channel considered. 

In Figs. [H 03 and [6] we show the scaled radial matrix 
elements for all the dipole allowed transitions between 
s, p, d and / Rydberg states in Cs and Rb. The radial 
matrix elements between high lying Rydberg levels are 
strongly peaked for states of similar energy so in practice 
only small positive or negative values of n s — m occur. 
The calculations reveal that in all cases the matrix ele- 
ments are large for at most three values of n s — rii. This is 
important since it limits the number of Forster channels 
which must be taken into account for an accurate calcu- 
lation of the interaction strength. The matrix elements 
are close to their asymptotic n 2 ao scaling for n > 50 and 
fine structure dependent effects are generally small ex- 
cept for the s — > p and p — > d transitions in Cs, which 
has a larger fine structure splitting than Rb. 

The energy defects S were calculated using rece ntly 
measured values for the Rb quantum defects [H, l39j |. 
and older data for Cs[I(|. Since the radial integrals scale 
as n 2 and the energy defects scale as n~ 3 the interaction 
strength usually scales as C§ ~ n 11 . As we will see below 
the asymptotic n 11 scaling is often broken for specific 
values of 71 < 100 where the quantum defects conspire 
to give near resonant Forster interactions. These special 
values of n may be particularly useful for engineering 
strong interactions without needing to access very high 
lying states. 
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FIG. 4: (color online) Radial matrix elements divided by n 2 
for transitions ns 1 / 2 — * n a p 3 / 2 (filled circles) and ns 1 / 2 —> 
n s Pi/2 (filled boxes) in Cs and Rb. 



1. ns-i/2 + ns 1 / 2 <-» n s pj + n t pj 

The first particular case is the excitation of ns\/2 states 
which can Forster couple to pairs of n s pj , n t Pj states with 
j = 1/2, 3/2. There are three possible fine structure chan- 
nels giving Forster defects 



81 (n s ,n t ) 



E(n sP3/2 ) + E(n t p 3/2 ) - 2£ , (ns 1/2 )(50a) 
E(n sP3/2 ) + E(n tPl/2 ) - 2E(ns 1/2 )(50b) 
E(n sPl/2 ) + E(n tPl/2 ) - 2 E(ns 1/2 ) (50c) 



The situation for Cs and Rb is shown in Fig. [7] 
The behavior is similar for both with the energy de- 
fect decreasing like 1/n 3 although the larger fine struc- 
ture splitting in Cs separates the pi/ 2 and p 3 / 2 chan- 
nels as compared to Rb. At n — 70 in Cs the 
strongest channels are (70,69), <5 2 (70,69), <5 2 (69,70), 
and <5 3 (70,69), giving C* 6 = 716,315,381,227 GHz /jm 6 . 
The next contribution is that from S%(71, 68)/2tt = 
—2.8 GHz which has much smaller radial matrix ele- 
ments giving Cq = 0.05 GHz ^m 6 . For Rb we find 
that 5i(70,69),<5 2 (70,69),(5 2 (69,70), and <5 3 (70,69) give 
C 6 = 799,543,589,437 GHz /1m 6 . The next contri- 
bution is that from <5i(71, 68)/2tt = -2.69 GHz which 
has much smaller radial matrix elements giving C§ = 
0.06 GHz /1m 6 . As we will see in Section HVl the vari- 
ation of the energy defects between channels results in 
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FIG. 5: (color online) Radial matrix elements divided by n 2 
for transitions npi/ 2 — * n s d 3 / 2 (open boxes) and nps/2 — * 
n sdz/2 (filled circles) in Cs and Rb. In Cs the additional 
fine structure transitions np^/2 — * n a d^/2 are also shown with 
filled circles and have slightly smaller matrix elements than 
the transitions to 6(3/2 for n s = n and n s = n — 1, and slightly 
larger matrix elements than the transitions to d 3 / 2 for n s = 



n — 2 and n. 

d 5 / 2 are less than 1% and are not shown 



n — 3. In Rb the differences between d 3 / 2 and 



an almost isotropic interaction for Rb, and some slight 
angular variation for Cs. 



2. npj + npj 



sSl/2 + n t s 1 /2 



The next case is coupling of pj states with j — 1/2, 3/2 
to S1/2 states. There are two fine structure channels with 
energy defects 

5t = E(n s s 1 / 2 ) + E(n t s 1/2 ) - 2E(np 3/2 ) 
§2 = E(n s s 1 / 2 ) + E(n t s 1/2 ) - 2E(np 1/2 ) 

which are shown in Fig. [5] 

For Cs with n s = n + 1, nt = n channel 1 has a reso- 
nance at n = 42 where Si/2ir = 15.7 MHz and the cor- 
responding interaction strength is Cq — —432 GHz /an 6 . 
The high n interaction strength in this channel is Cq — 
-2920 GHz /mi 6 at n = 70. The S 2 channel for n s = 



is substantially weaker giving C, 



-324 GHz /xm 6 at n 
a high n resonance for n s 



6 — 

70. The S 2 channel also has 
2, n t = n — 1 with 
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FIG. 6: (color online) Radial matrix elements divided by n 2 



for transitions nd 



3/2 



s/5/2 (open boxes) and nd 5 



/a 



(filled circles) in Cs and Rb. The additional fine struc- 
ture transitions nd 5 / 2 n sj 7/2 differ by less than 0.1% from 
the / s /2 case. 



S 2 /2tt = —5.77 MHz at n = 83. However the matrix el- 
ements are small so we get a relatively weak interaction 
of C 6 = 104 GHz /mi 6 . 

For Rb with n s = n + l,n t — n channel 1 has a res- 
onance at 11 = 38 where Si/2ir — —4.1 MHz the cor- 
responding interaction strength is C§ — 843 GHz /mi 6 . 
At n = 70, n s = n + 1, n t — n, the Si,S 2 channels give 
C 6 = -2820 and -767 GHz /mi 6 . 

We see that both species have resonances which pro- 
vide a strong interaction at relatively low n in addition 
to very strong interactions at high n. Unfortunately the 
strongest interaction occurs in the S\ channel which has 
angular zeroes (see Table A} so it is only useful for special 
geometries where the zeroes can be avoided. 



3. npj + npj <-» n a Si/2 + ntdj 2 

The next case is coupling of pj states with j = 1/2, 3/2 
to S1/2 and dj 2 states with j 2 — 3/2, 5/2. There are three 
fine structure channels with energy defects 

61 = E(n s s 1/2 ) + E(n t d 5/2 ) - 2E(np 3/2 ) 
S 2 = E(n s s l/2 ) + E(n t d 3 / 2 ) - 2E(np 3 / 2 ) 
83 = E(n s s l/2 ) + E(n t d 3/2 ) - 2E(np 1/2 ) 
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FIG. 7: (color online) Forster energy defects for nsi/2 
n s pj + ntPj coupling in Cs and Rb. 



FIG. 8: (color online) Forster energy defects for npj 
n s s i/2 + n tSi/2 coupling in Cs and Rb. 



which are shown in Fig. [9J The radial matrix elements 
are large for n s = n, n + 1 and rit = n — 2,n — 1, so we 
focus on these cases. 

For Cs the three cases n s = n, rit = n — 1, n s = n + 1, 
n t — n — 2 and n s = n + 1 , n t = n— 1 lie within a factor of 
10 in strength for all fine structure channels. At n = 70 
we find for Cq{5, n s , nt): 

C 6 (5x, 70, 69) = 111; C 6 (<5 2 , 70, 69) = 109 

C 6 (5 3 ,70,69) = 137; C 6 (Si, 71, 68) = 8.95 

C 6 (c5 2 ,71,68) = 8.2; C 6 (5 3 , 71, 68) = 11.0 

C 6 (5i, 71, 69) = -71.5; C 6 (S 2 , 71, 69) = -72.9 
C 6 (<5 3 ,71,69) = -59.0 

all in units of GHz yum 6 . 

In Rb the behavior is similar except the n s = n + 1, 
nt = n— 1 case has a larger energy defect than the others. 
At n = 70 we find for Cq(5, n Sy n t ) : 



C 6 (Sx, 70, 69) = 218; C 6 (5 2 , 

C6(6 3 ,70,69) = 253 Cfe($t, 

C 6 (5 2 ,71,68) = 61.0 C 6 (S 3: 

C 6 (5!,71,69) = -52.4 C 6 (5 2 , 
(76(53,71, 69) = -48.2 



70,69) = 217 
71,68) = 61.9 

71.68) = 71.0 

71.69) = -52.6 



channel does not and is therefore a good candidate for 
blockade experiments. 



4- npj + npj «-» n s dj 1 + ntdj 



The next case is coupling of pj states with j = 1/2, 3/2 
to dji states with j' = 3/2,5/2. There are four fine- 
structure channels with energy defects 

c5i = E(n s d 5/2 ) + E(n t d 5/2 ) - 2E(np 3/2 ) 

6 2 = E(n s d 5/2 ) + E(n t d 3/2 ) - 2E(np 3/2 ) 

S 3 = E(n s d 3/2 ) + E(n t d 3 / 2 ) - 2E(np 3/2 ) 

<5 4 = E(n s d 3/2 ) + E(n t d 3 / 2 ) - 2E(np 1/2 ). 

In Cs the strongest interactions occur for n s = n — 
1, rit =71—1, and n s = n — 2, n t = n — 1 and the 
corresponding energy defects are shown in Fig. |TD] At 
n = 70 we find for Ce(S,n s ,n t ) : 



(7 6 (5i,69,69) = -428. 
C 6 (6 3 ,69,69) = -478. 
<7 6 (5i,68,69) = 8.42; 
<7 6 (<5 3 ,68,69) = 7.75; 



C 6 ((5 2 ,69,69) = -451. 
C 6 (<5 4 ,69,69) = -334. 
C 6 (<5 2 ,68,69) = 8.36 
C7 6 ((5 4 ,68,69) = 10.3 



all in units of GHz ^m 6 . As can be seen in Table U the 
8i,5 2 channels suffer from zero eigenvalues, but the ($3 



all in units of GHz /im 6 . There is also a resonance at 
n = 68 for n s = n + l,n t = n — 3 for which the en- 
ergy defect is remarkably small, 8i/2ir = 0.57 MHz. This 
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FIG. 9: (color online) Forster energy defects for npj <-» 
risSi/2 + Tttdji coupling in Cs and Rb. For the cases with 
negative energy defect the defect is smallest in magnitude for 
the S3 channel, and for the the case with positive energy de- 
fect the 82 channel has the defect with smallest magnitude. 
The Si and 82 channels are very close together due to the 
small fine-structure splitting of the d states. 



value should be considered approximate since the quan- 
tum defects are not known accurately enough to predict 
the resonance to better than a few MHz. The small- 
est of the other resonances at the same value of n is 
62/271 = —250 MHz so to a good approximation we can 
consider just the first channel. The radial matrix ele- 
ments are R^ d — —352 ao and R^ d = —553 do giving 
C*6 = —63.5 GHz fim 6 . It should be noted that because 
the defect is so small the crossover to van der Waals be- 
havior occurs at R c = 8.7 /im despite the relatively weak 
strength of the interaction. 

In Rb the strongest interactions also occur for n s = 
n — l,n t = n — 1, and n s = n — 2, n t = n — 1 . At n = 70 
we find for Cq(S, n s ,n t ) : 

C 6 (<5i,69,69) = -112.; C 6 (5 2 , 69, 69) = -113. 

C e (5 3 , 69, 69) = -113.; C 6 (<5 4 , 69, 69) = -104. 

C 6 (<5i,68,69) = 48.6; C 6 (6 2l 68, 69) = 48.6 

C e (5 3 , 68, 69) = 48.0; C 6 (£ 4 , 68, 69) = 53.1 
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FIG. 10: (color online) Energy defects for npj+npj <-> n 3 dj 1 + 
ntdj 2 coupling in Cs and Rb. 



5- nd 3 /2, 5 /2 + ^3/2,5/2 <-» n-sPi/2,3/2 + ratPi/2,3/2 

The final groups of cases to consider are the d 3 / 2 , £^5/2 
states. Coupling of dj «-> occurs for four possible 
channels 

<5i = E(n s p 3/2 ) + E(n t p 3 / 2 ) - 2E(nd 5/2 ) 
82 = E(n s p 3/2 ) + E(n t p 3 / 2 ) - 2E(nd 3/2 ) 
S 3 = E(n s p 3/2 ) + E(n t pi/ 2 ) ~ 2E(nd 3/2 ) 
64 = E(n s px/ 2 ) + E(n t pi/ 2 ) - 2E(nd 3/2 ) 

which are shown in Fig. 111! 

In Cs the strongest channels are n s — n + l,n t = n + 1 
and n s = n + 2,n t = n + 1. At n = 70 we find for 
C 6 {S,n s ,n t ) : 

C 6 (Si, 71, 71) = 504.; C 6 (<5 2 , 71, 71) = 563. 

C e (8 3 , 71, 71) = 464.; C 6 (&i, 71, 71) = 394. 

C 6 (£i, 71, 72) = -11.1; C 6 (5 2 , 71, 72) = -10.2 

C 6 (S 3 , 71, 72) = -13.2; C 6 (<5 4 , 71, 72) = -13.5 



all in units of GHz /im 6 . 



all in units of GHz /im 6 . The same channels dominate 
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FIG. 11: (color online) Energy defects for nd 3 / 2 ,5/2 + 
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n sPi/ 2 ,3/2 + mpi/2,3/2 coupling in Cs and Rb. 
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FIG. 12: (color online) Energy defects for nd^/2,5/2 + 
^3/2,5/2 <-> «sPi/2,3/2 +ntfs/2,7/2 coupling in Cs and Rb. 



in Rb where we find 

C (5 lt 71, 71) = 132.; C 6 (S 2 , 71, 71) = 133. 

C7 6 (<5 3 , 71, 71) = 127.; C 6 (S 4 , 71, 71) = 122. 

C 6 (£i, 71, 72) = -65.4; C 6 (6 2 , 71, 72) = -64.8 

C 6 {5 3 , 71, 72) = -70.5; C 6 (6 4 , 71, 72) = -71.7 

in units of GHz ^m 6 . In this case there is a large dif- 
ference between the species with the total interaction 
strength summed over the four channels about eight 
times larger in Cs than Rb. However, all of these channels 
have zero eigenvalues and are therefore not immediately 
useful for blockade. 



6. nd 



'3/2,5/2 



+ WG(3/2,5/2 «-» ^sPl/2,3/2 + ^t/s/2,7/2 



The next case is coupling oindj +ndj <-> n s pj 1 +n t fj 2 . 
There are four possible channels 

6i = E(n s p 3/2 ) + E(n t f 7/2 )-2E{nd 5/2 ) (54a) 

5 2 = E(n s p 3/2 ) + E{n t f 5/2 )-2E(nd 5/2 ) (54b) 

5 3 = E{n s p 3/2 ) + E(n t f 5/2 ) - 2E{nd 3/2 ) (54c) 

5 4 = E{n sPl/2 )+E(n t f 5/2 )-2E{nd 3/2 ) (54d) 

which are shown in Fig. 1121 



-190. 
-192. 



In Cs the channels n s = n + l,n t = n — 3 and n s = 
n + 1, nt = n — 2 have comparable interaction strengths. 
At n = 70 we find for Cq(S, n s , n t ) : 

C 6 (£i,71,67) = 60.5; C 6 (<5 2 , 71, 67) = 60.5 

C 6 (5 3 ,71,67) = 64.3; C 6 (<5 4 , 71, 67) = 60.3 

C 6 (<Ji,71,68) = -190.; C 6 (6 2 , 71, 68) 

C 6 (£ 3 ,71,68) = -177.; C 6 (<5 4 , 71, 68) 

in units of GHz pm 6 . There are other cases such as n s = 
n + 4, n t = n — 5 which have resonances, e.g. 64/21: — 
—39 MHz at 71 = 65. However the radial matrix elements 
are too small to be useful. 

In Rb the strongest channel at large n is n s = n + 
l,nt = n — 1, with n s — n + 2, n t — n — 2 contributing 
about 30% as large a Cq and n s = n + 3,n t — n — 3 being 
substantially weaker. At n = 70 we find for Cq{8, n s ,nt) ■ 



C7 6 (<S 1 ,71,69) = 
C 6 (<5 3 ,71,69) = 
(78(^,72,68) = 
C 6 (<5 3 ,72,68) = 



-2530; C 6 (<5 2 ,71,69) = 

-2280; C 6 (<5 4 ,71,69) = 

-677; C 6 (<5 2 ,72,68) = 

-547; C 6 ((5 4 ,72,68) = 



-2530 
-3740 
-676 
-2330. 



in units of GHz ^m 6 . 

The n s = n + 2,n t = n — 2 channel is particularly in- 
teresting as it has a near resonance at n — 43 where 
Si/2tt = -8.3 MHz and 8 2 /2n = -6.0 MHz as well 
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FIG. 13: (color online) Energy defects for 71^3/2,5/2 + 
^3/2,5/2 «-» ^s/5/2,7/2 + ^i/5/2,7/2 coupling in Cs and Rb. 
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FIG. 14: (color online) Angular dependence of the resonance 
and blockade shifts for the 70s!/ 2 states in Rb and Cs at 
R = 9.2 fim. The C% parameters for Rb are given in the 
text and for Cs the values are: Cei ~ 712, Ce2 ~ 687, and 
6*63 = 213, in units of GHz /1m 6 . 



15.1, C 6 ((5i,67,68) = 188.0, and C 6 (5i, 68, 68) = -50.5, 
in units of GHz fan 6 . The other channels have similar 
strengths within about 15% of the given values. 

For Rb the strongest cases are n s — n t = n — 2, n s = 
n — 2,n t = Ti—l, and n a = n% = n— 1. At n — 70 we find 
for C 6 (S,n s ,n t ) : C 6 {8 U 68, 68) = 7.42, C 6 (5i, 68, 69) = 
77.7, and C 6 (5 1 , 69, 69) = -114., in units of GHz ^m 6 . 
The other channels have similar strengths within about 
1% of the given values. 



as n — 58 where 5^/271 = —6.9 MHz and n = 59 
where 64/211 = 8.6 MHz. The corresponding interaction 
strengths are C 6 (5i,45,41) = 391, C 6 (5 2 , 45, 41) = 539, 
C 6 (S 4 , 60, 56) = 6090, and C 6 (<5 4 ,61,57) = -5680, in 
GHz /1m 6 . The n = 58 resonance is the strongest we have 
found for any states with n < 70 and has a crossover 
length of R c — 11.8 /im. Unfortunately the interaction 
has angular zeroes (see Table H]) and is only useful for 
special geometries as discussed in Sec. IIVI 

7. nd 3 / 2 ,5/2 + nd 3 / 2 ,5/2 «-* n a f 5 / 2 j/2 + 714/5/2,7/2 

The final case is coupling of ndj ^ n s fj 1 +n t fj 2 . There 
are four possible channels to consider 

8 1 = E(nJ 7/2 )+E(n t f 7/2 )-2E(nd 5/2 ) (55a) 

h = E(n s f 7/2 ) + E{n t f 5/2 )-2E(nd 5/2 ) (55b) 

c5 3 = E(n s f 5/2 )+E(n t f 5/2 )~2E(nd 5/2 ) (55c) 

84 = E(nJ 5/2 )+E(n t f 5/2 )-2E(nd 3/2 ). (55d) 

The energy defects for these channels are all very similar 
due to the smallness of the fine-structure splitting of the 
d and / states. 

For Cs the strongest cases are n s — n t = n — S, n s — 
71 — 3, n« =n — 2, and n s — n t — n — 2 as shown in Fig. 
H31 At n = 70 we find for C 6 (8, n s , n t ) : C 6 (S U 67, 67) = 



IV. EFFECTIVE ANGULAR STRUCTURE OF 
THE FORSTER INTERACTION 

The choice of Rydberg states for blockade experiments 
is dictated by the strength and angular structure of the 
Forster interaction. A three dimensional distribution of 
atoms includes pairs with arbitrary relative orientations 
so that laser fields with laboratory fixed polarizations will 
generally couple to all possible two-atom eigenstates, in- 
cluding those with weak interactions. Excitation of these 
Forster zero states can be avoided either by choosing Ry- 
dberg states with near isotropic interactions, or by using 
carefully chosen interaction geometries. We give some 
representative examples of both approaches in this sec- 
tion. 

The prime example of a near isotropic interaction is 
the channel S\/ 2 + si/ 2 — > p + p which, as can be seen 
from Table [H is fully isotropic provided the fine structure 
of the p states is ignored. Accounting for fine structure 
gives the van der Waals Hamiltonian 

ffvdw = ^[C7 6 idiag(44,68,36,44) 

+ C 62 diag(28,4,36,28) 

+ C7 63 diag(8,32,0,8)]. (56) 

The Hamiltonian has been expressed in the basis (f|, 
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Tl+XT ^ LUJ ; where the arrows denote the electron 
spin projections for the two atoms and C6j is the coeffi- 
cient of channel j in Eqs. ([50|) . For the Rb 70si/2 state 
we find C 6 i = 794, C 62 = 1125, and C 63 = 42 7, and 
eigenvalues {891, 862, 862, 853} in units of GHz ^m 6 /i? 6 . 
Excitation of the Rydberg state 70s 1 / 2 |TT) in the labo- 
ratory frame gives the angular dependence shown in Fig. 
[T4l We see that even for Cs which has a relatively large 
fine structure splitting the interaction strength is close 
to isotropic. Because of this the resonance shift D and 
blockade shift B for each species are almost exactly the 
same. 

The angular average of the blockade shift in Rb can be 
approximated by the convenient expression 



B 70s = 1 MHz x 



9.77 ^m 
R 



(57) 
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FIG. 15: (color online) Angular dependence of the resonance 
and blockade shifts for the 43d 5/ / 2 state in Rb assuming a 
one- dimensional atom cloud with a — 3.0 /jm. 



which shows that a strong blockade is possible for two 
atoms with separations up to about 10 /^m. In appli- 
cations to ensembles containing atom pairs with a dis- 
tribution of R values we can calculate the blockade 
shift analytically using a Gaussian description for the 
atomic density. In a spherically symmetric cloud with 
radial density variance 3cr 2 the probability distribution 
is P(r) = (2T7a 2 )-^ 2 e- r2 / 2a2 . Replacing k£ w /I>£ by an 
angular mean (k 2 /D 2 } the spatially averaged blockade 
shift is 



The blockade shift then takes the form 



1 



N 



in 2 



K 

D 2 



(N 

N I k 2 \ (3.785 a) 



Hr dr' P{r)P(r') |r - r' 



r 2 



(58) 



where N is the number of atoms and we have used 
8648640 1 / 12 3.785. For the s 1/2 + s 1/2 -> p+p channel 
the interaction is isotropic and k 2 /D 2 = 9/16. Note that 
the strong weighting of the integrand towards large |r— r'| 
implies that the assumption of the van der Waals form 
for all molecular separations holds approximately even 
though at small separation the interactions may transi- 
tion into the resonant Forster regime. 

A second spatial distribution of interest is a quasi 
one-dimensional ensemble in an optical trap created by 
tightly focused laser beams. If T IQ \, the atomic temper- 
ature relative to the peak depth of the confining poten- 
tial, is small we can use a quadratic approximation to the 
potential about its minimum which leads to a Gaussian 
distributed density. To be specific consider a far off reso- 
nance trap (FORT) created by focusing a single Gaussian 
beam of wavelength A to a waist w (1/e 2 intensity radius) . 
When w is at least several times larger than A the trap 
provides a quasi one-dimensional distribution with prob- 
ability density P(z) — (2na 2 )~ 1 / 2 e~ z l 2a where a = 

whr^rel • FORT traps were used in several recent ex- 
periments with highly localized atomic clouds (20l. [111. |42| . 
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B 2 



N{n 2 /D 2 ) t 



(N- 



1)C 6 2 



dzdz'P(z)P(z'){z- z') 



t\12 



N 



(N-l) 



kT\ (3.0567 a) 12 
D 2 J e CI ■ 



(59) 



where (k 2 /D 2 ) s gives the interaction strength when the 
trap is tilted by an angle 6 with respect to the quantiza- 
tion axis z of the light, and we have written 665280 1 / 12 = 
3.0567. A single beam FORT with T re] = 2.5, w = 2.5 /zm, 
and A = 1.03 /im which is close to the parameters of our 
recent experiment (2Cj| gives a — 3.0 /im so 3.0567 a — 
9.2 \xm and the averaged interaction strength for 70s is 
that shown in Fig. [T3J 

In contrast to the s states the other channels in Table U 
exhibit strong angular effects. A case of particular inter- 
est is 43ds/2+ 43d 5 / 2 — * 45p3/ 2 +41fj, which is within a 
few MHz of being Forster resonant for Rb. The Forster 
defects are -6.0, -8.3 MHz for j — 5/2,7/2, which give 
Cq coefficients of C%\ = 391, 6*62 = 539 GHz^tm 6 with 
Cgj referring to channel j in Eqs. ([54")) . The effective 
Cq, averaging over the degeneracy of the two channels, 
is C 6 = |C 6 i + |C 62 = 454 GHz^m 6 . If we were to 
ignore angular effects we would naively expect a strong 
interaction of order 454/860 = 0.53 times that shown in 
Fig. [TJ] but at a much smaller value of n which relaxes 
the laser power requirements for fast excitation. 

However, as mentioned above, the 36 D v coefficients 
for the d 5 1 2 + d b / 2 — * P3/2 + / channel cover a huge range, 
with two M = states being extremely small. To be 
explicit assume we start in the 87 Rb / = 2, nif = hy- 
perfine ground state and use 7r-polarized excitation light 
which couples to the Rydberg state (expressed in a basis 
aligned with the light polarization) 
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FIG. 16: (color online) Angular dependence of the resonance 
and blockade shifts for the A3d 5 / 2 state in Rb for excitation 
of the | f , f ) Rydberg state and assuming a one-dimensional 



atom cloud with a = 3.0 fim. 



where the second two numbers in each ket are the mag- 
netic quantum numbers for the nuclear spin. The latter 
are conserved in the excited state (assuming the hyper- 
fine interaction can be neglected there). We must there- 
fore average 1/B 2 for each of the four terms to calculate 
the angular dependence shown in Fig. [15] The small 
value of B for angles between the light polarization and 
the molecular axis near 90 deg. render this interaction 
a poor choice for blockade in a spherical ensemble. We 
see that even at the optimum angle the blockade shift is 
only 65 kHz which is about 20 times smaller than in Fig. 
[TH Furthermore, the resonance shift D is larger than 
the blockade shift, and has a somewhat different angular 
structure. 

It is worth noting that it is possible to choose interac- 
tion geometries which largely avoid the small D v coeffi- 
cients of 0.003, 0.002 for this channel. The state with the 
smallest D v is interchange symmetric (the eigenvector 
can be found in the accompanying EPAPS material) 
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where the kets are given in the form \mk,mi). The state 
with the second smallest D v is very similar to the above 
but is interchange antisymmetric. Coupling to \iPf) can 
be strongly suppressed by using a + excitation light with a 
one-dimensional trap aligned along z, and Zeeman select- 
ing a target m,j level to excite the two-atom state ||, |). 
The angular distribution for this state is shown in Fig. 
[TBI We see that for a trap aligned parallel to the quan- 
tization axis the blockade strength is 0.25 MHz, which 
is more than 8 times larger than when exciting the state 
given in Eq. (|60[) . Finally, we note that a similar trick can 



be used to render the 58d 



3/2 ■ 



58d 



3/2 



60pi /2 + 56/ 5/2 

resonance discussed above usable for blockade in one- 
dimensional geometries. The blockade shift obtained at 



9 = by exciting ||, |) is a very large \B\ = 2.9 MHz in 
a trap with a = 3 /im. 



V. CONCLUSIONS 

In this paper we have considered in detail the effects 
of Zeeman degeneracy on blockade experiments relying 
on van der Waals interactions to allow only single-atom 
excitations. The figure of merit for blockade is sensitive 
primarily to the weakest interactions between the various 
degenerate Rydberg states. For many convenient Ryd- 
berg states, the degeneracies result in particular linear 
combinations of Zeeman sublevels having zero or nearly 
zero van der Waals interactions. This problem can some- 
times be avoided using special geometries and choices of 
light polarization, but care must be taken. 

We have catalogued the long-range potential curves for 
a large number of angular momentum channels likely to 
be of interest to blockade experiments, with sufficient in- 
formation to allow researchers to quantitatively evaluate 
the van der Waals interactions for a wide range of exper- 
imental situations. 

All blockade experiments reported to date have used 
samples whose spatial extents are substantially larger 
than the range of the van der Waals interactions. In these 
situations, dipole-dipole or van der Waals interactions 
can play a dominant role in suppressing Rydberg exci- 
tation under conditions where quantum blockade would 
not be possible. 
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APPENDIX A: 



QUADRUPOLE-QUADRUPOLE 
INTERACTION 



There is also a quadrupole-quadrupole (Q-Q) interac- 
tion that contributes a R~ 5 term in the long-range po- 
tential Hi: 



Vqq 



R 5 



{Qa <S> Qb} 40 = 



'70 



P 



(Al) 

where the atomic quadrupole moment operator is Q p = 
er 2 y 'Air /5Y2 P . As with the van der Waals interaction, 
the quadrupole-quadrupole interaction causes transitions 
between different Zeeman levels. 

The following argument shows that for most cases the 
quadrupole-quadrupole interaction will be considerably 
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smaller than the van der Waals. For the Q-Q interaction 
to dominate over the van der Waals interaction, we need 



e 2 (r 2 r e 4 (r) z 
> 



(A2) 



R 5 5R e 

The matrix element factors are roughly equal, so we find 

„2 



R > 



350 /im 



(A3) 



for a 1 GHz (or smaller) value of S that is usual for Ry- 
dberg states. Thus the Q-Q interaction should be negli- 
gible. 



APPENDIX B: DIPOLE-DIPOLE INTERACTION 
IN THE COUPLED BASIS 

Instead of the basis | ja'^aib'mb) we could equally well 
do the calculations in the coupled basis \j a jbJM). Us- 
ing multipole tensors and recoupling algebra, a simple 
relation for the matrix elements of the dipole-dipole in- 
teraction can be obtained. 

The multipole tensors are defined using the tensor 
product formalism of 32] (3.1.7): 

T*jL = ® {J\}.JM = E C™JM- m \jrn){j'M - m| 

rn 

(Bl) 



where the time reversed bra is defined as 



{jm\ = (-iy +m (jfn\ 



(B2) 



The spherical component of the position operator r of an 
electron can, for example, be written in terms of T\ : 



3s J 



(B3) 



as can be verified by taking matrix elements of both sides 
of the equation. 

The dipole-dipole interaction is proportional to the 
spherical tensor 



Us\\r\\j){jt\\r\\j) 



r} 
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x{{|i s )®{i|}i®{|i t )®{i|}i} 20 (B4) 

where we are assuming that in the initial states the elec- 
trons have the same angular momentum j, and the cou- 
pling of interest is isolated to a single state where the 
electrons on atoms a and b have angular momenta j s and 
jt- 

We can now use recoupling algebra to rewrite this in 
terms of the coupled states {\j s ) £g> \itL} KM = \j s jtKM) 
and {\j) ® \j)} JM = \jjJM), using |1 (3.3.2 (11)): 



{a (g) b} 2 



(js\\r\\j)(jt\\r\\j) 



}{\K)®{J\} 



\jsjl 

KJ [ KJ2 . 
(jsjl 

[KJ2 



KJ 



where [J] = \/2J+ 1. 

In this basis, it is often found that the van der Waals 
eigenstates are heavily weighted with a single value of 
J, suggesting that in many cases J is an approximately 
good quantum number. 
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